SBW/UNp02/Figures V.l 

Resonant X-ray Scattering Study of Magnetic and 
Electric-quadrupole Order in Uo.75Npo.25O2 

S. B. Wilkins/'2 J. A. Paixao,^^ R. Caciuffo/'^ P. Javorsky/ F. Wasting 
J. Rebizant/ C. Detlefs,^ N. Bernhoeft,^ P.Santini,^ and G. H. Lander^ 

^European Commission, Joint Research Centre, 
Institute for Transuranium Elements, 
Postfach 2340, Karlsruhe, D-76125 Germany 
^European Synchrotron Radiation Facility, 
Boite Postal 220, F-38043 Crenoble CEDEX, France 
^Physics Department, University of Coimbra, Coimbra, 3004-516 Portugal 
^ Dipartimento di Fisica ed Ingegneria dei Materiali e del Territorio, 

Universitd Politecnica delle Marche, 1-60131 Ancona, Italy 
^Department de la Recherche Fondamentale sur la Matiere Condensee, 
CEA, F-38054 Grenoble CEDEX, France 
^ Dipartimento di Fisica, Universitd di Parma, 
Viale Parco delle Scienze, 1-43100, Parma, Italy 
(Dated: February 2, 2008) 



Abstract 

We have used element specific X-ray resonant scattering to investigate the M edge resonances 
in a single crystal of Uo.75Npo.25O2. Earlier neutron diffraction and Mossbauer studies had shown 
the presence of long-range AF order below Tjv = 19 K, with sizeable magnetic moment both on 
the U and the Np ions. RXS results confirm the presence of an ordered dipole magnetic moment 
on the Np ions, but also reveal the presence of an anti-ferro arrangement of the electric quadrupole 
moments on both U and Np ions, with the same propagation vector that defines the magnetic 
structure. Prom the azimuthal dependence of the intensities we are able to determine the exact 
configuration of the quadrupolar ordering. The intensities of Bragg peaks associated with magnetic 
dipole and electric quadrupole order have different temperature dependences. On cooling below 
Tjv, the magnetic dipole order develops faster on the uranium ions, with magnetic order on the 
Np ions increasing at a lower rate. At the same temperature, quadrupolar order on both the U 
and Np ions occurs along with an internal Jahn-Teller lattice distortion. 

PACS numbers: 75.25.+Z, 75.10.-b 
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I. INTRODUCTION 



Actinide dioxides, An02, have fascinating physical properties whose understanding is 
still a challenge, despite half a century of extensive investigations carried out by both the- 
oreticians and experimentalists. The complexity of their physical behavior stems from the 
subtle interplay between crystal field (CF), superexchange interactions and electron-phonon 
coupling, which results in peculiar multi-polar ordered-phase diagrams. 

At room temperature An02 systems have the fee crystal structure of CaF2, with space 
group Fm3m and tetravalent An'^^ ions. Ordered magnetism below Tjv = 30.8 K was sug- 
gested in UO2 by the presence of a large anomaly in the temperature dependance of the 
heat capacity^, and confirmed by magnetic susceptibility^ and neutron diffraction^ mea- 
surements. Neutron spectroscopy was soon after applied to elucidate collective excitations, 
vibrational and magnetic, in the low temperature phase*". Evidence of a strong spin-lattice 
interaction prompted a theory introducing the concepts of Jahn- Teller (J-T) effects and elec- 
tric quadrupole interactions in the actinide materials^. Some years later, neutron diffraction 
actually revealed an internal J-T distortion of the oxygen sub-lattice in UOa^, and eventually 
it was demonstrated that the magnetic structure is 3-k type, i.e. Fourier components of all 
three members of the star of = (001) are simultaneously present on each magnetic site of 
the lattice.-^ These Fourier components of the magnetic moment, /x^, are perpendicular to 
the propagation vector k, so that the structure is transverse. The overall symmetry remains 
cubic, and the magnetic moments /xq point along the [111] directions of the cubic unit cell. 
As Hk can assume two mutually orthogonal directions for each k value (for instance, Hk par- 
allel to [010] or [001] for k = [100]), there two distinct 3-k transverse magnetic structures 
(Figure 1). 

Experiments on powder samples with chopper spectrometers at spallation neutron sources 
allowed a precise and reliable determination of the CF Hamiltonian acting on the U ions^, 
and experiments with triple-axis spectrometers and axial polarization analysis have been 
performed to study collective excitations in the ordered phase, and dispersive magnetic 
excitons in the paramagnetic phased. The experimental picture is therefore exhaustive, but 
a satisfying theoretical description of low-temperature collective excitations is still missing. 

Magnetic order was expected to occur also in Np02, for which a large A-like anomaly in the 
heat capacity^^ and a maximum in the magnetic susceptibility curved were observed at Tq = 
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FIG. 1: a) Shows a 3d representation of the 3-fc magnetic structure. Uranium atoms are shown as 
blue spheres with the corresponding dipole magnetic moment direction shown as an arrow. Oxygen 
atoms are shown as red spheres, b) The longitudional structure as projected onto the a — b 
plane, i.e. on moving from an atom at (000) with a moment along [111] to the atom at (^^0) the 
moment direction changes to [TTl]. c) and d) are transverse structures. Following the reasoning 
above the atom situated at (^^0) has moment directions [Til] and [111] for c) and d) respectively. 



25 K. However, neither Mossbauer spectroscopy^^ nor neutron diffraction experiment o^^i^^i^^ 
fiave found evidence of magnetic ordering in tliis compound. Tlie experiments set an upper 



4 



limit on the Np ordered magnetic moment of 0.01 hb, ^ value that is surprisingly small given 
the Kramers nature of Np^"*" ions (5/^, ^1^/2 in Russell-Saunders approximation), and the 
large effective Curie- Weiss paramagnetic moment deduced from susceptibility measurements 

Indeed, Resonant X-ray Scattering (RXS) experimentsii' at the Np M4 edge exclude 
usual magnetic dipole ordering at Tq, but provide direct evidence of long-range order of an 
electric quadrupole moment with Fs symmetry. The phase transition is purely electronic 
and involves neither internal nor external crystallographic distortions, so the symmetry of 
the system remains cubio^^*^^^. It has been suggested that the primary order parameter is 
associated with T\ magnetic octupoles, ordering in a 3-fc longitudinal structure with (001) 
propagation vector. The ground state (GS) is a singlet with no static dipole magnetic 
moment in which magnetic octupole order breaks time reversal invariance and induces an 
electric quadrupole as the secondary order paramete r^°i^^ . 

The magnetic behavior of Ui_xNpx02 solid solutions has also been extensively inves- 
tigated by powder neutron diffraction, ^^'''Np Mossbauer spectroscopy, and magnetic sus- 
ceptibility measurement o^^ i^^i^"^ . The ^^^Np Mossbauer isomer shift indicates that Np ions 
are tetravalent, with 5/^, ^1^/2 electronic configuration, independently of temperature and 
composition. Magnetic order has been detected in the composition range from x = 0.15 
and X = 0.75, with transition temperatures decreasing with increasing x, from T/v ~ 17 K 
(x = 0.25) to T/v ~ 11 K (a; = 0.75), showing that substitution of Np for U reduces the 
strength of U-U exchange interactions. For x = 0.25, the magnetic structure is transverse 
antiferromagnetic type I, as in UO2, with propagation vector (001), and an average ordered 
magnetic moment fi = 1.8{1)^b- As Mossbauer data are compatible with an ordered mag- 
netic moment fi^p = 0.6(2)/iB on the Np ions, this gives fiu = 2.2(2)/!^ on the U ions. For 
larger Np content, the magnetic order is short range and incommensurate^^. 

II. RESONANT X-RAY SCATTERING 

Resonant X-ray Scattering (RXS) occurs when a photon is virtually absorbed by exciting 
a core electron to empty states, and subsequently re-emitted when the excited electron and 
the core hole recombin o^^i^^ . This process introduces anisotropic contributions to the X-ray 
susceptibility tensoi^, the amplitude of which increases dramatically as the photon energy 
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is tuned to an atomic absorption edge. In the presence of long-range order of magnetic 
moments, electronic orbitals occupancy, or spatially anisotropic valence-electron clouds, the 
interference of the anomalous anisotropic scattering amplitudes may lead to the excitation 
of Bragg peaks at positions forbidden by the crystallographic space group. For absorbing 
atoms belonging to the actinide series, both L- and M-edges are of practical use, as their 
energy is large enough for Bragg diffraction to be allowed in most actinide compounds. The 
experiments reported here have been performed at the U and Np M4 and M5 edges, involving 
respectively 8^3/2 5/ and 8(^5/2 ^ 5/ electric-dipole. El, virtual excitations. 

The amplitude of magnetic scattering from a given atom at the 3d — 5/ threshold is 
the product of a geometrical factor, depending on the polarization of the incident (e) and 
scattered (e') photons, and a linear combination of resonance strength factors, F^^m? 

FWoc(6xe')[Fi,i-Fi,_i] (1) 

Magnetic scattering therefore requires a rotation of the photon polarization, and a differ- 
ence in the transition probabilities to intermediate states with M = ± 1. This asymmetry 
can arise because of the net spin polarization of the 5/ states, or from a difference between 
overlap integrals, resonant energy, or lifetime for spin-up and spin-down channela^&iSSiSS. 

RXS can also signal long-range order of anisotropic charge distributions. It is the as- 
phericity of the atomic electron density that generates the anomalous tensor component in 
the atomic scattering factor. If the electron clouds on different sublattices have different 
orientations, superlattice reflections occur due to the reduced translational symmetry. This 
is the so-called Templeton or ATS (Anisotropic Tensor Susceptibility) scattering^S and is 
observed, for instance, in case of antiferro-electric order of the quadrupole moments. The 
scattering amplitude arising from El transitions can be described by second-rank tensors 
/(r), invariant under the point symmetry of the scattering atom. If f{Q) is the Fourier 
transform of f{r), the scattering amplitude can be written as 

FPl(Q)(xe'-/(Q)-e, (2) 

where e and e ' are the polarization vectors of the incident and scattered beam, where the 
rotation about the azimuth is implicitly included in the electric field vector. 
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III. EXPERIMENTAL DETAILS 



Uo.75Npo.25O2 single crystals have been grown using the chemical-transport reaction 
method at the Institute for Transuranium Elements in Karlsruhe. Starting materials were 
obtained by stoichiometric co-precipitation by ammoniac, calcination under air at 1073 K, 
followed by reduction under Ar/H2 (5%) at 1073 K. Pellets were fabricated and sintered at 
1973 K for 6 hours in Ar/H2 atmosphere. Charges of about 5 g of starting material were 
crushed and welded in a quartz tube with TeCU (4.5 mg/cm^) as transport agent in a 10~^ 
torr vacuum. The feed substances were maintained at 1323 K and transported under the 
influence of the temperature gradient into the growth region at 1223 K for 14 days. The sin- 
gle crystals obtained were characterized by chemical analysis, scanning electron microscopy, 
single crystal x-ray diffraction and x-ray diffraction on a powder obtained from crushed 
crystals. No differences in the average composition or evidence for clustering effects have 
been observed by these techniques. However, differences in the spatial homogeneity of the 
Np ions distribution in individual crystals grown by the method described above cannot be 
excluded. 

Two specimens with a mass of 5.20 mg and 3.63 mg have been selected to measure the 
heat capacity as a function of temperature and magnetic field amplitude, using a PPMS-9 
Quantum Design platform. The RXS experiment has been performed on the 1D20 beamline 
of the European Synchrotron Radiation Facility (ESRF) in Grenoble, France. For these 
measurements we used a crystal 0.204 mm^ in volume, with a mass of 2.26 mg and the 
[111] axis perpendicular to the exposed face. The sample was of excellent quality with a 
mosaic width of 0.05 °. For safety requirements, the crystal was encapsulated in a sealed 
double-wall container with beryllium windows. 

The ID20 spectrometer exploits the fundamental harmonic emission of two 42-mm-period 
undulators. A double crystal Si(lll) curved monochromator provides sagittal focusing and 
selects the appropriate X-ray energy. Vertical focusing and filtering of higher harmonics is 
achieved by Si mirrors, that deliver at the sample position a 99 % linearly polarized beam, 
with polarisation perpendicular to the scattering plane (a) and incident wavevector kj. A 
mosaic crystal positioned between sample and detector was used to analyze whether the 
polarization of the scattered beam (with wavevector k/) is parallel (tt) or perpendicular (cr) 
to the vertical scattering plane. 
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A closed-cycle refrigerator equipped with an azimuthal-rotation stage provided a base 
temperature of 10 K and the possibility to perform azimuthal scans. The spectrometer was 
operated in the photon energy interval from 3.5 to 4.0 keV, including the U and Np M- 
edges (U: Em„ = 3.552 keV, Em^ = 3.728 keV; Np: Em, = 3.664 keV, Em^ = 3.845 keV), 
using the (111) reflection from a Au analyzer (Bragg angle of 43.2° at the Np M4 edge, 
and 41.6° at the U M4 edge; no correction was performed to account for the departure 
from the ideal 45° orientation of the analyzer crystal). The energy resolution was about 
0.5 eV (full width at half maximum), and the beam flux at 3.845 keV was estimated as 
5 X lO^^photons s~^ . Integrated intensities Im for Bragg reflections were obtained numerically, 
by fitting a linear background and a Lorentzian squared lineshape to the measured peaks. 
Corrections for the size of the footprint of the incident beam and for absorption were applied 
as Jc = /m /^(l + sin a/ sin/3), where a and (3 are the angles formed by the incident and 
scattered beam directions with the sample surface, and /i is the energy dependent absorption 
coefficient. Because of the strong absorption all reflections must be found in Bragg geometry. 



IV. RESULTS OF THE EXPERIMENT 



Figure 2 shows the heat capacity as a function of temperature measured for the 
Uo.75Npo.25O2 sample with 5.20 mg mass. Data were collected from 0.7 to 300 K in zero 
fleld and under an external magnetic fleld of 9 T. The base temperature was limited by 
self-heating effects, that become sizeable below 3 K. Contributions to the measured heat ca- 
pacity from the sample platform (including grease used to establish thermal contacts) were 
measured separately and subtracted from the total signal. The results reveal the presence 
of a A-like cusp at 18.9 K, and two small anomalies at 16.8 K, and 15 K. The application of 
a 9 T magnetic fleld does not produce appreciable effects; the maximum at 18 K undergoes 
a very small shift to a lower temperature while its amplitude slightly decreases. The speciflc 
heat of the isostructural compound Th02 is assumed to be a close estimate of the vibrational 
contribution, and is shown as a solid line in Figure 2. As previously reported for both UO2 
and Np02ii, the magnetic contribution is appreciably different from zero at temperatures 
far above the anomaly also in the solid solution. 

The broken line in Figure 2 represents the results of measurements performed on a second 
single- crystal with the same nominal composition of the flrst one. The curves obtained for 
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the two samples are quite different; for the second one, the main anomaly occurs at 16.8 K, 
whereas small anomalies are visible at 18.4 K and 15 K. The characterization analysis 
suggests the same stoichiometry within experimental errors, and X-ray diffraction does not 
show line broadening, thus excluding major clustering effects. On the other hand, the subtle 
interplay between magnetic and electric multipolar interactions in this system could be very 
sensitive to even small changes in the chemical short-range order. The large difference 
observed in the heat capacity of the two samples could therefore reflect small variations in 
the probability distribution of Np neighbors around U ions. 

The presence in the diffraction pattern of resonant superlattice reflections with propa- 
gation vector (001), forbidden in the paramagnetic phase, proves that long-range order of 
some multipolar moment characterizes the low-temperature ground state. The integrated 
intensity of the (112) Bragg peak at 10 K is shown as an example in Figure 3, as a func- 
tion of the photon energy across the U and the Np M5 and M4 absorption edges. These 
data were collected without analyzing the polarization of the diffracted beam, and therefore 
represent the sum of the intensities in the a ^ a and a ^ vr polarization channels. The 
spectrum can be reproduced as the superposition of four oscillators that interfere with each 
other. The solid line in Figure 3 is obtained by assigning to each oscillator a Lorentzian 
lineshape centered at the El electric-dipole energy, including corrections for self-absorption 
and energy dependence of the beam attenuation. 

At 10 K, a resonant behavior has been observed for both the o" ^ vr and the a ^ cr 
components of the spectrum. This is shown in Figures 4 and 5, around the U and the Np 
M4 edge, respectively. Note that the a ^ a have Lorentzian squared lineshapes, and are 
shifted by —2 eV wrt. the cr — vr peak. This has been previously observed in systems 
where the a — > a signal is known to originate from the F'^l term.^^ The ratio between the 
intensities of the a a and the cr — vr peak is about 30% for Np, and about 3% for 
U. For strictly magnetic dipole ordering, there is no signal in the a ^ a channel. Thus, 
the finite ratio indicates that the structurally forbidden peaks do not stem uniquely from 
dipole magnetic-order. Our RXS results are compatible with dipole magnetic-order together 
with an ordered distribution of electric quadrupole moments. Resonant scattering from 
orientationally ordered aspherical electron clouds appears in both polarization channels. On 
the other hand, the shape of the resonance across the M4 absorption edge is different for the 
a ^ a and the cr ^ tt contributions. Whilst the a ^ vr peak can be fitted to a Lorentzian 
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FIG. 2: Heat capacity measurements performed with a Quantum Design PPMS-9 system. Data 
were taken in an applied field of (circles) T and (triangles) 9 T. A A-like peak and two anomalies 
are clearly visible. The solid line is a fit to heat capacity data measured for the non-magnetic iso- 
structural compound Th02. The broken line corresponds to results obtained for a second sample 
with the same nominal composition. 

function, the a — > a resonance has a Lorentzian-squared hneshape. A Lorentzian peak is 
usually observed when the intensity enhancement is due to dipole magnetic order, whereas a 
Lorenztian-squared peak has been observed in Np02, both in the rotated and the unrotated 
polarization channeU^. 

The narrow energy profile of the Bragg intensity at the Np M4 edge in Np02 has been 
associated with a core-level structure, and attributed to magnetic and charge con- 

tributions to the El Bragg amplitude described by Np 5f multipoles of rank 3 and 42^. 
This explanation cannot be applied in the present U + ions in UO2 do not carry 
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FIG. 3: Normalized integrated intensity as a function of energy for the (112) reflection of 
Uo.75Npo.2502at 10 K. Polarization analysis was not performed. The solid line is a simulation 
assuming four dipole oscillators, each represented by a Lorentzian function. 

magnetic octupoles. On the other hand, in ref. the RXS signal observed in Np02 is 
interpreted entirely as a consequence of electric-quadrupole order, both for the a ^ a and 
the a —>■ 7T channels, as the magnetic ordered moment in Np02 is vanishingly small. The 
results presented here seems to support these latter interpretations. 

As the magnetic ordered moment in Np02 is vanishingly small, in that case the reso- 
nant scattering is entirely due to electric-quadrupole order. It therefore appears that the 
Lorentzian squared shape is a characteristics of quadrupole order in these compounds. We 
conclude that in this solid solution, the a ^ n signal is dominated by the magnetic dipole 
contribution, whilst the a —>■ a resonance is due only to the electric-quadrupole order. 

The results of \1/ (azimuthal) scans about the (112) magnetic peak, taken at T = 10 K with 
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FIG. 4: Integrated intensity of the (112) superlattice reflection, at 10 K, as a function of photon 
energy in the vicinity of the U M4 absorption edge. Upper panel o" ^ vr lower panel a ^ a. The 
solid line (upper panel) is a fit to a Lorentzian lineshape with a FWHM of 8 eV. In the lower panel 
the solid line represents a fit to a Lorentzian squared lineshape with the same width and a sloping 
background. 

photons energy tuned to the Np and U M4 edges, are shown in Figure 6. The intensity of the 
peak is measured while the sample is rotated about the scattering vector Q, kept constant at 
the chosen value. The origin of the azimuthal angle, is defined where the [lll]-direction 
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FIG. 5: Integrated intensity of the (112) superlattice reflection, at 10 K, as a function of photon 
energy in the vicinity of the Np M4 absorption edge. Upper panel cr ^ vr lower panel a ^ a. 
The solid line (lower panel) show a fit to a Lorentzian squared lineshape. The a ^ tt contribution 
is dominated by the magnetic dipole signal and exhibits a Lorentzian lineshape (upper panel). 
Intensities are normalized to the maximum of the a ^ vr signal at the M4 dipole threshold. Note 
that the cr ^ vr signal (upper pannel) is strongly influenced by the dominant signal from the U M4 
resonance, see Figure 
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is within the scattering plane. Oscillations with two-fold symmetry are observed for both 
the U and Np signals in the a a channel. On the other hand, the a — vr intensities are 
independent of \E'. The azimuthal dependence of the peak intensity is similar at the U and 
Np M4 edges, indicating the same electric quadrupole orientation on the two ions. 

The solid line in Figure IHl is a simulation of the scattered intensity for a 3-A; anti- 
ferro quadrupole structure assuming a incoherent addition of the two transverse s-domains 
(See Figure P). The lower panel of Figure IHl shows a simulation of the longitudional 3-fe 
anti-ferroquadrupole structure. This model clearly does not fit the data. Details of the 
calculations used to obtain the simulated azimuthal dependence are given in the appendix. 

Figure 7 shows the temperature dependence of the (112) integrated intensities, normalised 
to unity at 10 K. The reported data were taken at the M4 dipole threshold of either Np or 
U, for both polarization channels. Intensity was collected as a function of time while the 
sample temperature was varied at a rate of 2 K min~^. The cr — vr intensity is dominated by 
the magnetic dipole contribution and is therefore proportional to the square of the ordered 
magnetic moment. The a ^ a channel intensity, on the other hand, is proportional to the 
square of the electric quadrupole order parameter. 

On cooling, the magnetic-dipole order develops first on U at T/v ~ 19 K; the order- 
parameter increases with an apparent plateau between T/v and 18 K. The magnetic order on 
Np, and electric quadrupole order of both ions, is estabhshed at T/v but the corresponding 
order-parameters increase at a lower rate. Note that the AF phase transition in UO2 is 
strongly discontinuous, whereas in the solid solution the quadrupolar and the magnetic 
phase transition are continuous. The temperature dependence of the U a ^ vr and a a 
intensity is remarkably different, whilst for Np the corresponding difference is marginal. 

Besides resonant superlattice reflections, diffraction peaks corresponding to a deformation 
of the oxygen cube have also been found. The angular positions of the observed peaks 
corresponds to those expected for the inhomogeneous internal distortion observed in UO^^, 
and resulting from combinations of normal modes of the oxygen sub-lattice. For these 
measurements, we used a photon energy of 7.5 keV and a Ge (333) analyser to select the 
a ^ a polarisation channel, where structural Bragg peaks give a contribution. As an 
example, we show in Figure 7 (solid line) the intensities measured for the internal lattice 
modulation measured at the (421) reflection. It appears that the onset of the Jahn- Teller 
distortion coincides with the onset of the a ^ a scattered intensity, i.e. the long range order 
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FIG. 6: Azimuthal dependence of the integrated intensity of the (112) superlattice reflection in the 
a ^ a channel at both the Np M4 (top panel) and U M4 (middle panel) absorption edges. Lines 
are intensities calculated assuming a electric qi^gdrupole moments along the (111) direction, for a 
addition of two noncollinear transverse antiferro-quadrupolar domains (see Figure^. The lower 
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FIG. 7: Scattered intensity of the forbidden reflections as a function of temperature. The data for 
the (112) reflection were taken at the M4 edges of U and Np in both polarization channels: (filled 
circles) U o" — > vr intensity; (open circles) \J a a intensity; (upper triangles) Np cr — vr intensity; 
(down triangles) Np a — s- a intensity. The solid line intensities are measured for the (421) internal 
lattice modulation at 7.5 keV, indicating the onset of a Jahn- Teller distortion at Tjy. Solid lines 
are fit to a power law. 



of the electric quadrupolar moments. 

A departure from the cubic symmetry was also checked by following angular position and 
width of the (006) lattice reflection, but no evidence for a lowering of symmetry was found. 
This confirms the 3-k nature of the ordering. On the other hand, as shown in Figure 8, an 
anomalous volume contraction of the cubic cell is observed on cooling between T/v and 18 K. 
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FIG. 8: Temperature dependence of the lattice parameter of the C/o.75-^Po.2502 soHd solution. Error 
bars are smaller than the symbols. A anomalous contraction can be observed in the temperature 
interval between 18 and 20 K. 

V. DISCUSSION 



The CF spectrum of actinide dioxides is known with high degree of precision from inelastic 
neutron scattering (INS) experiments. The CF ground state of U in UO2 is a Fs triplet, 
whereas the CF ground state of Np in Np02 is a Fg quartet. The CF potential changes little 
among UO2, Np02, and even PUO2. Therefore, the CF acting upon the U and Np ions in 
solid solutions UxNpi_x02 is expected to be the same as in the pure compounds. 

The F5 triplet carries magnetic-dipole and electric-quadrupole degrees of freedom. U-U 
dipolar and quadrupolar interactions are provided by superexchange (SE) and additional 
U-U quadrupolar interactions are carried by phonons through the so-called virtual-phonon 
exchange. In UO2 dipolar interactions are stronger than quadrupolar interactions and would 
produce on their own a second-order AF transition at around 25 K. A quadrupole moment 
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is induced by dipoles as secondary order parameter (OP). In turn, this quadrupole induces 
a lattice distortion of the oxygen cage. Although quadrupolar interactions are too weak 
to produce a purely quadrupolar phase transition, they contribute to the free energy of 
the system when quadrupoles are induced by dipoles. This contribution is believed to be 
responsible for the change of the magnetic transition from second- to first-order with an 
increase of Tjv up to about 31 K. The AF structure, type-I transverse 3-A;, is probably 
stabilized by quadrupolar interactions. As yet no theoretical model has been able to explain 
the occurrence of the 3-k structure. 

In Np02, the Fg quartet carries magnetic-dipole, electric-quadrupole and magnetic- 
octupole degrees of freedom. Np-Np dipolar, quadrupolar and octupolar interactions are 
provided by superexchange (SE) and additional Np-Np quadrupolar interactions are carried 
by phonons. The latter are expected to be appreciably smaller than the corresponding U-U 
interactions because of the smaller quadrupolar polarizability of Np ions. Dipole-dipole in- 
teractions are expected to be smaller as well?^. Octupolar interactions appear to dominate, 
and induce a second-order AF transition with no dipole moment, but with an associated 
quadrupole secondary OP. No lattice distortion occurs^^. 

Ion-ion interactions in solid solutions UxNpi_x02 are particularly rich. As a reason- 
able guess, U-U and Np-Np interactions may be assumed to be similar to those in pure 
compounds. The lack of a scaling framework for SE makes U-Np interactions difficult to 
estimate a priori. A few qualitative observations are formulated in the following from a 
virtual-crystal-approximation- like perspective, i.e. focusing on "average" properties and ne- 
glecting effects specifically associated with the disorder in the sample. 

Since the magnetic phase transition of Uo.75Npo.25O2 occurs at a lower T/v than in UO2, Np 
acts as a magnetic diluent, and therefore U-Np dipolar interactions must be much smaller 
than the corresponding U-U interactions. Since adding Np makes the transition second- 
order, quadrupolar interactions must be smaller than in UO2. U-Np octupolar interactions 
are absent since octupoles are quenched in U. Thus, the phase transition of Uo.75Npo.25O2 
is driven mainly by U-U dipolar interactions. U quadrupoles (and the associated lattice 
distortion) follow U dipoles as secondary OPs but their contribution to the free energy is 
too small to affect the phase transition qualitatively (i.e. to change the value of T/v and the 
order of the transition). Their slow rate of growth below T/v reflects a large (3 exponent, 
which is typical of a secondary OP linearly coupled to the square of the primary OP (in 
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mean-field (5 is twice the primary's OP j3). 

Np ions are expected to play a more passive role. The leading interaction of Np02, i.e. 
octupole-octupole, has no significant effects in Uo.75Npo.25O2 since it involves the minority 
Np-Np pairs but not U-Np and U-U pairs. As U ions order, Np ions feel the developing U 
dipoles and quadrupoles through U-Np interactions, and themselves acquire ordered dipoles 
and quadrupoles. The delay by which these Np moments are detected by X-rays can be 
understood quite naturally within this picture: since the U-Np dipolar interaction by which 
the U dipolar order is transmitted to Np ions is (relatively) weak, thermal fluctuations on 
Np ions are quenched by the phase transition more slowly than on U ions. So, Np moments 
appear at T/v together with U moments, but they grow more slowly as T decreases, although 
eventually as T — > the full moments must be recovered. 

In order to check this qualitatively, we performed calculations for the solid solution based 
on the virtual-crystal-approximation (VCA)2i. Russell-Saunders coupling was assumed, and 
for simplicity only the CF ground multiplets of U and Np were included in the calculation, 
neglecting therefore excited CF states. Although the wavefunctions of the Fg quartet of 
Np depends on the cubic CF parameter x, results do not depend qualitatively on x (x = 
—0.75 was used^°). The VCA calculation assumed a 3-fc structure and contained four order 
parameters to be determined self-consistently at each temperature, i.e. the size of the 
dipole moment along (111) and of the axial quadrupole moment along (111) for U and Np, 
^111(^^)5 ^111 (Np), (5iii(U), Qiii(Np). There are six parameters, i.e. the dipole-dipole 
(J) and quadrupole-quadrupole {K) interaction strengths for U-U pairs, for U-Np pairs 
and for Np-Np pairs. Np-Np couplings play a minor role since these pairs are in minority. 
The two leading dipole-dipole couplings J(U-U) and J(U-Np) are determined by the known 
transition temperatures of pure UO2 and of the solid solution. If quadrupole-quadrupole 
interactions are assumed to be mainly lattice-transmitted, then i^(U-Np) and i^"(Np-Np) 
can be obtained by scaling K(U-U) by using the known quadrupolar polarizability of Np 
ions. K{\]-\]) is the parameter most difficult to determine, because the role of quadrupole 
interaction in the structure of pure UO2 is not quantitatively clear. 

As an example, we show in Figure IHl calculated RXS intensities for a parameter set with 
a fairly small value of -ft"(U-U), i.e. about half the value which had been tentatively used 
in UOs^. {3 exponents are mean-field like in the VCA, (3 = 0.5 for U and Np dipoles and 
(3 = 1 for U and Np quadrupoles. The curves are in qualitative agreement with experiment, 
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FIG. 9: Calculated temperature dependence of the normalized RXS intensity from U and Np 
dipoles and quadrupoles for Uo.75Npo.25O2. 

with U dipoles ordering more quickly than all other moments. This remains true if K{U-U) 
is varied. Increasing K{\J-\J) mainly increases the curvature of the red squares in the figure, 
i.e. it accelerates the ordering of U quadrupoles. However, U dipoles always prevail near 
T/v, because of their smaller P exponent. 

As far as the specific heat C is concerned, a quantitative understanding of the measured 
curve is difficult for at least two reasons. The first is that the sample dependence of C sug- 
gests that slight inhomogeneities in the sample may have pronounced effects on C. Besides 
being difficult to quantify and characterize, inhomogeneities cannot be included in a VGA 
calculation, as disorder is averaged out in the VGA. The second reason is that the measured 
C implies that the magnetic entropy S at T/v is remarkably reduced with respect to what 
expected {S = A;b(0.25 ln(4) + 0.75 ln(3)) per actinide ion), while the measured 5* is lower 
than ks ln(2) per actinide ion (entropy reduction is observed in pure UO2 and Np02 as well, 
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although the reduction is smaller). Thus, it is a priori impossible for any mean- field model 
such as the VGA to reproduce the measured C quantitatively. A quahtative feature of C 
whose origin is understandable within the VGA is the presence of slight oscillations below 
T^r. These seem to refiect the onset of contributions to the free energy from the "delayed" 
ordering of Np ions. These contribute little to the main peak of C at Tjv, and produce 
features below T/v when their ordering begins to become appreciable. 

The VGA provides a qualitatively satisfactory theoretical framework for Uo.75Npo.25O2, 
yet there is one point (besides the entropy reduction problem) which is difficult to under- 
stand. This is the low value of the saturation ordered moment of Np ions. There appears 
to be no way to obtain a moment as low as 0.6/iB with the observed 3-fc structure.— This 
problem would remain even if theoretical approaches different from the VGA were used, 
as it is associated with the structure of the single-ion GF states of Np. These appear to 
carry a moment along (111) always much larger than 0.6/iB, no matter whether the GF is 
varied or whether quadrupole interactions are included. In particular, since the dipole and 
quadrupole order parameters appearing in the 3-fc structure commute, the bare GF value of 
the moment cannot be reduced by quadrupolar interactions. This would not be true if the 
structure was 1-k. 

The experiments on the mixed oxides have shown that at the 25% Np doping (x = 0.25) 
the transverse nature of the 3 — transition predominates, despite the fact that in pure 
Np02 the ordering of the quadrupoles is longitudinal 3 — fc. Although the agreement in 
the top frame of Fig. 6 is not outstanding, it must be remembered that these experiments 
were performed on an off-specular reflection [the (112) from a (111) face] and this creates 
more difficulties experimentally than when a specular reflection is examined, as was the 
case in pure Np02.— The crucial aspect is that the maximum and minimum for azimuthal 
scans agree with the calculations for a transverse 3 — A; quadrupolar ordering and completely 
disagree for the longitudinal model. This change across the solid solutions, from transverse 
in UO2 to longitudinal in Np02 also gives a possible mechanism for the unusual effects 
found in the range 0.50 <x < 0.75 [Refs. 1330 as with increasing Np concentration 
increasing frustration can be expected between these competing quadrupolar interactions. 
It will therefore be particularly interesting to try and observe the quadrupolar ordering in 
the compositions with higher Np content. 

Finally, with these observations of the quadrupolar ordering in the mixed oxides we note 



21 



in each case the Lorentzian squared shape of the resonances and their displacement to lower 
energy by about 2 eV (see Figs. 4 and 5). They arise from the term in the RXS 
cross section2^i2& and have also been observed in the work we reported on high-order terms 
in multi— fc structures.— The strong similarity between these signals and those found^^ in 
pure Np02 suggests that in Np02 as well the principle RXS signal arises from a quadrupole 
signal, and not a more complicated process, as suggested recently. 

Acknowledgments 

The authors wish to thank Matt Longfield for his help in the early part of these exper- 
iments. SBW and PJ would like to thank the European Commission for support in the 
frame of the 'Training and Mobility of Researchers' programme. Part of this work was made 
possible thanks to the support of the European Community-Access to Research Infrastruc- 
tures action of the Improving Human Potential Programme (IHP) in allowing access to the 
Actinide User Laboratory at the Institute for Transuranium-Karlsruhe under the contract 
HPRICT-200100118. 

APPENDIX: AZIMUTHAL DEPENDENCE 

In a cubic 3-k structure the complete star of A; = (A;00) is simultaneously present in each 
volume element. Therefore there is only one k-domain. However, several orientations of the 
Fourier components of the magnetic moment relative to the corresponding wave vector are 
possible. In particular, there are one longitudinal and two transverse directions. 

The longitudinal orientation is unique, so for this structure there is only one S-domain. 
This would be the structure of Np02 if there was a magnetic dipole moment. In the trans- 
verse case, there are two possible orientations which may be converted into each other by 
applying crystallographic symmetries. This leads to 2 equivalent but distinct S-domains. 
This is the case in UO2. 

Below, we calculate the quadrupolar term of the cross section) scattering from these 
3 constructions. 

We assume the crystal structure is FCC with 4 atoms per unit cell. For each atom in turn 
we can calculate the moment direction at this position for the longitudinal and transverse 
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Atom Long. Trans. Domain A Trans. Domain B 

(0,0,0) /^i = (1,1,1) /^i = (1,1,1) /ii = (1,1,1) 

(iiO) /X2 = (1,1,1) /^2 = (1,1,1) /li = (T,l,T) 

(iO,i) /X3 = (T,1,T) /^3 = (1,1,1) /X3 = (1,1,1) 

(0,i,i) /X4 = (1,T,T) /X4 = (1,1,1) /14 = (1,1,1) 



TABLE L Moment directions on the four atoms of the unit cell for the single longitudinal and two 
transverse S-domains 

domains (See Figure P). The resulting directions are given in Table L 

To form the scattering tensor for quadrupole scattering we first evaluate the quadrupole 
operator Qij which in cartesian coordinates is given by 



Qij — l^i^j — -^ij y^Xf^kfJ'k) 

k 



(A.l) 



This matrix is symmetric and traceless. For all of the above structures /i^ = l,Vi, which 
leads to the general form 



Q 



(A.2) 



Evaluating this for the 3 S-domains from the moment directions given in Table I gives 
scattering tensors for each atom in the unit cell given by Table II for the quadrupole moment. 
From these the structure factor can be created as 



n 

which for the Bragg reflection (1, 1, 2) gives 

/(112) = Ti +T2 — T3 — T4 



(A.3) 



(A.4) 



yielding scattering tensors of the form 



(112) 



1 



1 
\ J 



(A.5) 
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Atom 



Long. 



Trans. Domain A 



Trans. Domain B 



(0,0,0) 



1 1 

1 1 
1 1 



1 1 

1 1 
1 1 



1 1 

1 1 
1 1 



i" 


1 




1 





T 




T 







T 




1 





1 




1 





1 1 
1 T 



(io,i) 



1 1 
1 T 



10 1 
T T 



1 1 
T 1 



1 1 
T 1 



10 1 
1 T 



1 1 
T T 



TABLE IL Scattering tensors for the four atoms in the FCC unit cell for the longitudinal and 
transverse cases. 



7(112) 



'^O 0^ 

1 

Vo 1 oy 

'^O 1^ 



1 

This can be evaluated to give a total structure factor of 



(112) 



V 



7 



p[2] 
^(112) 



e' ■ /(112) ■ e, 



(A.6) 



where e and e' are the polarization vectors of the incident and scattered beams, respectively. 
Here we consider only cr — > cr scattering, as the cr — > tt channel is dominated by magnetic 
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scattering (-F^) term. 

To evaluate this, first we construct a coordinate system of tfie crystal. We therefore, 
construct 3 unit vectors based on the reciprocal space vectors q (the scattering vector) and 
^azref (the azimuth reference vector), ci, £2 and C3. These are related by, 

= c, X (A.7) 

\ I "'azrci I / 
Ci = C2 X C3. 



We then do a rotation around £3 to transform into the coordinate system of Blume and 
Gibbs.— We therefore obtain 3 new unit vectors of the form, 

= C3 

U2 = Ci sin + C2 cos (A. 8) 

Ui = Ci COS (j) — C.2 sin 0. 

which again are linear combinations of the reciprocal space vectors g, fcazref and the cross 
product of these two. We can then define the polarization vectors for the incident and exit 
beams in the new coordinate system. For a polarization this is simply equal to the unit 
vector ■U2, i-e. the polarization vector is normal to the scattering plane (given by the cross 
product of the scattering vector and the azimuthal reference vector). For vr polarization 
the Bragg angle of the reflection has to be taken into consideration. This is included here 
for completeness but is not used for these calculations. We therefore have the polarization 
vectors in the rotated coordinate system given by 

= e'„ = U2 (A.9) 
= uisoiOb — U'^cosOb (A. 10) 

= — iii sin^^ — ^3 cos^^ (A. 11) 



Given these vectors equation IA.6I can be evaluated to give the amplitude of the scattering 
for each required value. 
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